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Abstract 

We study the scenario where a sneutrino is the next lightest supersymmetric par- 
ticle (NLSP) and decays into a gravitino and standard model particles. The daughter 
particles such as neutrinos and quarks cause electro- and hadronic showers and affect 
big-bang nucleosynthesis. It is found that despite a small branching ratio, four-body 
processes including quarks in the final state give the most stringent constraint on the 
sneutrino abundance for gravitino mass of 1 — 100 GeV. Pion production by high en- 
ergy neutrinos is important when the sneutrinos decay at ~ 1 sec. We also discuss the 
thermal leptogenesis in the sneutrino NLSP scenario. 



1 Introduction 



Among various possibilities, supersymmetry (SUSY) is a prominent candidate for physics 
beyond the standard model since it possibly solves various serious problems in particle physics 
and cosmology. In particular, from cosmology point of view, supersymmetric models contain 
a good candidate for dark matter; with i?-parity conservation, the lightest superparticle 
(LSP) becomes stable and can be dark matter. It is widely known that, if parameters are 
properly chosen, relic density of the LSP agrees with the dark matter density suggested by 



where h is the Hubble constant in units of 100 km/sec/Mpc. This fact, as well as other 
motivations of supersymmetry, provides a strong motivation to consider supersymmetry as 
a new physics beyond the standard model. 

If we consider supersymmetric models, however, several problems may also arise. In 
particular, in (local) supersymmetric models, superpartner of the graviton, i.e., gravitino, 
exists. Gravitino is a very weakly interacting particle and it may cause serious cosmological 
problems [2]. In this paper, we concentrate on the case where the gravitino is the LSP.*^ In 
this case, gravitino becomes a potential candidate for dark matter. 

If gravitino is the LSP, the next to the lightest superparticle (NLSP), which is assumed to 
be the lightest superparticle in the minimal supersymmetric standard model (MSSM) sector, 
decays into gravitino with very long lifetime.*^ Importantly, the lifetime may be longer 
than 1 sec. Thus, in the evolution of the universe, relic NLSPs may decay during and/or 
after the big-bang nucleosynthesis (BBN). With the decay of the NLSP, electromagnetic 
and hadronic showers may be induced. Energetic particles produced in the showers induce 
photo- and hadro-dissociation and production processes of light elements (D, ^He, ''He, 
^Li, and ^Li). Such processes may significantly change the prediction of standard BBN 
scenario and, consequently, resultant abundances of light elements may significantly conflict 
with observations PU]- Recent detailed analysis have shown that the model is seriously 
constrained when the NLSP is neutralino or stau [TH I12j. 

However, the lightest neutralino and stau are not the only possible candidates for the 
lightest superparticle in the MSSM sector. In particular, if the scalar neutrino becomes 
lightest among the superparticles in the MSSM sector, we expect that the constraints from 
BBN are drastically relaxed. This is because the dominant decay mode of the parent particle 
(i.e., sneutrino) is i> — ipp, and hence the daughter particles are both weakly interacting.*^ 

#^For recent study for the cases with unstable gravitino, see, for example, (See, also, ^01^1 [3 El El-) 
*Hn this paper, we assume i?-parity conservation. 

^^Onc may also consider the case where the sneutrino is the LSP and gravitino is the next LSP. Then, 
the primordial gravitino decays into neutrino-sneutrino pair. In this case, however, gravitino is likely to 
decay also into charged particles since there should exist charged slepton which is almost degenerate with 
the sneutrino. Thus, the most stringent constraints from BBN is expected to be from photo-dissociation 
processes induced by emitted charged particles; constraints on such cases can be read off from old studies. 
(See, for example, [ElIZllHl-) Thus, in this paper, we concentrate on the case where the gravitino is the LSP. 
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In order to discuss the BBN reaction in this class of set up, it is crucial to understand the 
effects of high-energy-neutrino injection on BBN. Such a study was partly done in |T2]; in [T^ . 
however, only photo-dissociation processes induced by high-energy neutrino injection was 
discussed, and hadro-dissociation processes were not considered. In addition, the sneutrino 
NSLP was studied in pjjj but the processes induced by the high energy neutrinos were 
ignored. 

In this paper, we consider cosmological constraints on models where (i) the scalar neutrino 
becomes lightest among the superparticles in the MSSM sector, and (ii) gravitino is the 
LSP. We pay particular attention to constraints from BBN. Organization of this paper is as 
follows. In the next section, we discuss properties of sneutrino which plays important role 
in our scenario. In Section El we discuss effects of long-lived sneutrino on BBN. Then, in 
Section HI we show results of our numerical study. Implication of our results for thermal 
leptogenesis scenario is discussed in Section Section IHl is devoted for conclusions and 
discussion. 

2 Properties of Sneutrino 

Before discussing cosmology, we first introduce the model we consider. In particular, we 
discuss properties of sneutrino which plays important role in our study. 

Here, we adopt the particle content of the MSSM. Important assumptions in our study 
are that the lightest superparticle in the MSSM sector is sneutrino 9, which is in SU{2)l 
doublet slepton L and that the lightest superparticle is gravitino ip. 

We first discuss when sneutrino can become the lightest superparticle in the MSSM 
sector. For this purpose, for simplicity, we adopt SU{5) grand unified model, which is 
strongly motivated by the fact that, with the particle content of MSSM, three gauge coupling 
constants meet with a good precision at the grand-unified-theory (GUT) scale Mgut — 
2 X 10^^ GeV. Then, left-handed (s)leptons and right-handed down-type (s)quarks are in 5 
representations of SU{5) while other (s)fermions are in 10 representations. 

Detailed mass spectrum of superparticles depends on the mechanism of SUSY breaking.*'^ 
In the case of gauge-mediated SUSY breaking jTHI, mass spectrum of MSSM superparticles 
are so constrained as far as the messenger sector respects SU{5). Consequently, sneutrino 
does not become the LSP in the MSSM sector. 

If we consider supergravity-mediated SUSY breaking scenario, on the contrary, soft SUSY 
breaking masses for sfermions arise from (non-renormalizable) interactions in Kahler poten- 
tial. Consequently, in this case, there exist large number of free parameters in the soft SUSY 
breaking terms. Thus, hereafter, we concentrate on supergravity-mediated SUSY breaking 
scenario. 

"^^In anomaly mediated SUSY breaking scenario |14| , gravitino becomes much heavier than superparticles 
in the MSSM sector and hence it cannot be the LSP. Thus, in anomaly-mediated models, the mass spectrum 
which we consider cannot be realized and hence our following analysis does not have relevance in such cases. 
Thus, in the following, we will not consider anomaly-mediated models. 
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When the soft SUSY breaking terms arise from Kahler interaction, soft SUSY breaking 
masses for 5 and 10 can be independent. In this case, all the sfermions in 10 representation 
can become heavier than sneutrino by pushing up the SUSY breaking masses for 10. In 
addition, it is usually the case that right-handed down-type squarks becomes heavier than 
left-handed slepton because of the renormalization group effects due to gluino mass, even 
though their SUSY breaking masses are degenerate at the GUT scale. Thus, it can be easily 
realized that left-handed slepton becomes the lightest among the sfermions in supergravity- 
mediated SUSY breaking scenarios. In "conventional" scenarios, however, it is usually the 
case that the bino mass becomes smaller than the sneutrino mass. This is due to renormal- 
ization group effects. Denoting the unified gaugino mass as M1/2 and SUSY breaking mass 
squared of 5 representation as m|, electroweak-scale value of the mass squared of SU{2)i 
doublet slepton is estimated as TB]*^ 

m|(M„eak) ^ mi + 0.53M2/2. (2.1) 

(Here, M1/2 and m| are defined at the GUT scale.) In addition, bino mass is given by 

Mi(M„eak) ^0.41Mi/2. (2.2) 

Thus, with the assumption of m| > 0, Ml becomes smaller than m|(Mweak) and sneutrino 
becomes heavier than bino. 

However, we can easily avoid this conclusion by assuming m| < 0. As we mentioned, 
some part of SUSY breaking parameters originate from direct interaction between SUSY 
breaking fields and observable-sector fields. Denoting the superfields for SUSY breaking 
field and observable-sector field as z and 0, respectively, the following interaction may exist: 

Ant = ^ j d'ez*z^*4>, (2.3) 

where A is a coupling constant. After SUSY is broken, scalar component in acquires 
SUSY breaking mass squared as — -r4-|-^^P, where is the F-component of z. Assuming 

iWpj 

A > 0, negative contribution to mass squared of scalar field is obtained. Importantly, A is a 
free parameter in Kahler potential and its sign is unknown from low-energy effective theory 
point of view. Thus, we assume that the A parameter for doublet slepton is positive and 
that m|(M„eak) < Of course, in this case, there exists true minimum where the doublet 
slepton acquires non-vanishing vacuum expectation value. This is not a serious problem 
since the transition rate to the true minimum is small enough so that the transition to the 
true minimum does not occur in the cosmic time scale. 

When the relation m|(Mweak) < ^1 holds, sneutrino may become the lightest super- 
particle in the MSSM sector in large fraction of the parameter space. Indeed, neglecting 
left-right mixing, masses of sneutrino u and left-handed charged lepton are given by 

1 

ml = m|(M„eak) + 2"^|cos2/?, (2.4) 
'^^Here, we have neglected the effect of Yukawa coupling constants; it may slightly decrease m|(Mwoak)- 
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m|(M„cak) + (sin^6'w 




- j m| cos 



(2.5) 



where is the Weinberg angle, and tan (3 is the ratio of vacuum expectation values of two 
Higgs bosons. Thus, when tan/? > 1, 9 becomes lighter than II. Thus, hereafter, we consider 
the case where the sneutrino is the lightest superparticle in the MSSM sector. In addition, 
as we mentioned, we consider the case where the LSP is the gravitino ip. In this class of set 
up, BBN may be affected by late-time decay of primordial sneutrino which freezes out from 
the thermal bath when the cosmic temperature drops below the sneutrino mass. 

Next, let us consider decay processes of z/. Decay of the sneutrino is dominated by 
the two-body process i) —>■ ipu; as we will see, decay rates for three- and four-body decay 
processes are much smaller than Tj^^^iy. Thus, the decay rate of sneutrino is given by 



where ^3/2 and mj> are the masses of gravitino and sneutrino, respectively, while ~ 
2.4 X 10^^ GeV is the reduced Planck scale. The lifetime of sneutrino t^, is calculated by 
using the two body decay rate given in Eq. ()2.6|) . In Fig. Q we plot contours of constant 
on 1713/2 vs. plane. Importantly, as the gravitino mass becomes smaller, the decay rate 
of j) becomes larger since the interaction with the longitudinal component of the gravitino is 
more enhanced. Thus, if the gravitino mass becomes smaller than ~ 0.1 — 1 GeV, lifetime of 
sneutrino becomes shorter than ~ 1 sec. In this case, primordial sneutrinos in early universe 
decay before BBN starts and no constraints are obtained from BBN. On the contrary, with 
larger gravitino mass, lifetime becomes longer and we should carefully consider hadro- and 
photo-dissociation processes induced by the decay of sneutrino. 

Even though the dominant decay of i) is via two-body process, it is also necessary to 
consider three- and four-body decay processes in the study of non-standard BBN reactions. 
In particular, processes like i) i:uZ^*^ and z> ^ ijlW^*\ followed by Z^*^ ^ // and 
_> ff'^ are important. (Here, the superscript "*" is for vertual particles. In addition, 
/ represents quarks and leptons, and / denotes charged lepton.) This is because charged 
and colored particles are produced only via three- and four-body decay processes. Those 
charged and colored particles efficiently interact with background particles and induce photo- 
and hadro-dissociation processes. As we will see in the following sections, these three- and 
four-body processes may significantly change predictions of standard BBN scenario in some 
region of parameter space. 

In order to study effects of decay modes with colored and/or charge particles in the 
final state, we calculate the decays rate for the processes u 'ip^ff, — ^ i^Uf, and 
u ipwy. Here, we approximated that the neutralinos and charginos are purely gaugino 
(or Higgsino). In addition, we neglect effects of Yukawa coupling constants. Then, for the 
processes i) — >■ ipvff and i) — > iplff, we take account of the Feynman diagrams shown in Fig. 
121 (In our analysis, we do not use the narrow-width approximation for the productions of "on- 
shell" Z- and IV-bosons; effects of "on-shell" Z and W productions are taken into account 
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Figure 1: Contours of constant lifetime of the sneutrino on vs. mi, plane. We have 
shaded the region with iriy < m^j2-, which we are not interested in. 

at the phase-space region where the invariant mass of // (or //') system becomes close to 
the gauge-boson mass.) In Fig. 13 we plot the contours of constant hadronic branching ratio, 
which is defined as 



■ 9 9 



(2.7) 



As one can see, the hadronic branching ratio is very small. As we will see, however, hadronic 
decay processes play important role in discussing dissociation of light elements. We also 
calculate decay rate for the process z/ — > ipv^., which occurs with the diagram like Fig. Wic) 
(without attaching // final state). We found, however, that the decay rate for this process 
is negligibly small with the parameter set we use in our numerical analysis. 

In order to study effects of photo-dissociation processes, we also calculate averaged "vis- 
ible energy" emitted from the decay of one sneutrino: 



1 



{u- 
vis 



f 



vis 



■'•Plfl') 



(2. 



where {E^'^^ ■*) are averaged energy carried away by charged particles and photon in in- 
dividual decay modes !>—>■■■. We have checked that the quantity {eI^^ ■*) are typically 
0(10 %) of the sneutrino mass when ^ n^'u- However, ii^vis is much smaller than rrii 

since the branching ratios for events with charged particle (s) in the final state are small. 
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v|/vf f\|/vf fv\|/f fV|/f fv 
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Figure 2: Feynman diagrams for the decay processes i) — >■ il'i'ff ((a) — (d)) and D — >■ 
■^^/Z' ((e) — (h)). Here, 5, and represent bine, neutral wino, and charged wine, 
respectively. 



3 BBN with Long-Lived Sneutrino 

Now, we are at the position to discuss effects of long-lived sneutrino on BBN. As we men- 
tioned, there are two types of decay modes of sneutrino; one is two-body decay mode while 
the other are three- and/or four-body ones. Reactions caused by these decay modes are 
different, so we discuss effects of these decay modes separately. 



3.1 Two-body decay 

First, we discuss effects of the dominant decay process u — > ipu. The gravitino produced in 
the decay is a very weakly interacting particle, so it is irrelevant for BBN. 

Neutrino, on the contrary, may affect abundances of light elements. Once emitted, the 
energetic neutrinos may scatter off background leptons via weak interaction. Consequently, 
several kinds of particles may be pair-produced. 

First, charged leptons may be produced via the following processes: 

i^ii^i,BG e"e+, (3.1) 
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Figure 3: Contours of constant hadronic branching ratio on m3/2 vs. plane. 



z^ii^i.BG (3.2) 
Vi,Ve,^G /W~e"^, (3.3) 
I^e^'/i,BG e~/i+, (3.4) 

where i = e, /i, r is flavor index, and the subscript "BG" is for background particles. The 
muons emitted in the above processes quickly decay into electrons and neutrinos. Thus, 
the above processes produce energetic electrons (and positrons) which cause electromagnetic 
cascade and energetic photons in the cascade induce photo-dissociation processes of light 
elements. Effects of these processes have been already studied in jT^]; it was pointed out 
that, in some part of the parameter space, abundances of light elements are significantly 
affected if energetic neutrinos are injected. 

Other possible effect is due to the production of pion pair which was not considered in ^3] ■ 
High energy neutrinos scatter off the background neutrinos and electrons and produce pions 
as 

ViVi^^G 7r"7r+, (3.5) 
i^jCi^BG 7r°7r^. (3.6) 

The nucleus-pion interaction rate is ~ 10* sec~^ x(T/MeV)^ which is larger than the decay 
rate of the charged pion (~ 4 x 10^ sec~^) for 1 MeV. Therefore, the charged pions 
produced at T ~ 1 MeV scatter off the background nuclei and change protons (neutrons) 
into neutrons (protons) via 

Tx^'p rnx^ or 77-7, (3.7) 
TT^n pn^ or p'j, (3.8) 
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which increases the n-p ratio and synthesizes more ^He. On the other hand, because of very 
short hfetime, the neutral pions decay before they scatter off the background nuclei. 

In order to estimate effects of the high energy neutrino induced processes, we have to 
numerically solve the Boltzmann equation describing the time evolution of the high energy 
neutrino spectrum taking account of neutrino-neutrino (charge lepton) scattering as well as 
production of charged leptons and pions. (The details of the Boltzmann equation and its 
solution will be presented elsewhere [H].) 



3.2 Three- and four-body processes 

Even though the branching ratio for the three- and four-body decay processes are much 
smaller than 1, such decay processes are very important since colored and/or charged par- 
ticles are directly emitted from these decay processes. Energetic colored and charged par- 
ticles may significantly change the prediction of standard BEN scenario. Effects of these 
colored and charged particles are classified into three categories: photo-dissociations, hadro- 
dissociations, and p ^ n conversion. 

Photo-dissociation processes are induced by energetic photons in electromagnetic shower 
which is caused by charged particles emitted from D. With given background temperature, 
the distribution function of energetic photons depends on total amount of energy injected 
by particles with electromagnetic interaction, and is insensitive to the energy spectrum of 
primary particles. Thus, once i^vis is obtained, the energy distribution of energetic photons in 
the electromagnetic shower can be obtained. Then, photo-dissociation rates are obtained by 
convoluting energy distribution function and cross sections of photo-dissociation reactions. 
For details of our treatment of photo-dissociation processes, see |.2Qj . 

For the study of hadro-dissociation processes, it is necessary to obtain energy distributions 
of (primary) hadrons which are produced after the hadronization of quarks emitted from 
u. As a first step to derive distribution functions of hadrons, we calculate invariant- mass 
distribution of qq^'^ system produced by z> — > ipuqq and z/ iplqq'. For example, for the 
neutral-current event u "ipi^qq, we numerically estimate the following quantity for each 
quark flavor q: 

dVi,_^yqq Nc f 2 7^ 7^ lA^ |2 



where A^'c = 3 is the color factor, Aif^^^^qq is the matrix element, and rriqq and m^qq are invari- 
ant masses of qq and uqq system, respectively. In addition, d^x,y represents (infinitesimal) 
two-body phase space of x and y. For each quark flavor q, spectra of hadrons are calculated 
by PYTHIA code [2I| . We denote the spectrum of nucleus A^ in the center-of-mass frame of 
qq, which is obtained after the hadronization of one pair of qq, as 



JN y-^N 1 ™qq) 



dN 



N 



dE^ 



(cm) 



(3.10) 
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where denotes the number of nucleus A^. It should be noted that /^™^ is not the spectrum 
in the rest frame of u. In order to take account of the effect of boost, we numerically calculate 
the averaged boost factor of the qq system for as a function of mqq: 

mqq j dmf,^^d^^^(^uqq)d^^^(^qq)d^q^q\Mu^^^qg\ 

where Egg is the energy of qq system in the rest frame of u. With this quantity, we estimate 
the spectrum of the nucleus N in the rest frame of 9, approximating that the boost factor 
for a fixed value of mqq is universally jqq{mqq). In addition, for simplicity, we approximate 
that the distribution of initial q (and hence q) jet is isotropic in the center-of-mass frame of 
qq. Then, we obtain 



dN, 



N 



dE 



N 



/ dcoseqdml^^^fi^-\E^-^-mqq)^^j^, (3.12) 



,^^.qq 2r,^^^qqJ ^ dE^ ' ^ ' dm^q 



where, in the above formula, -E^^™^ is related to En as 



En = %q [E^r^ + ^,,V^J,^"^'-^^cos^^,) , (3.13) 



with Pqq = y 1 — 7gg^. Effects of charged-current events z> — > iplqq' are also treated in the 
same way. 

We have calculated the spectra of p and n taking account of all possible neutral- and 
charged-current events. These hadrons cause hadronic shower and induce hadro-dissociation 
processes. In our analysis, in addition, we have also estimated the number of charged pions 
produced by the decay of i> with the same procedure as the case of p and n. Such charged 
pions become the source oi p ^ n conversion process, which changes the number of ^He. 
Once the spectra of hadrons are obtained, effects of hadro-dissociation and p ^ n conversion 
are studied following [3 |H] . 



4 Numerical Results 

Now, we are at the position to show our numerical results. In our analysis, we have followed 
the evolutions of the number densities of light elements by numerically solving Boltzmann 
equations. For this purpose, we have modified Kawano code ^22j, including non-standard 
processes discussed in the previous section: 

• p ^ n conversion by pions produced via vvbg ~^ tt+tt^ and three- and four-body 
decays of v. 

• Photo-dissociation processes by charged leptons produced via w-qg ~^ ■ 
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Hadro- dissociation processes by hadrons produced via three- and four-body decay pro- 
cesses of p. 



• Photo-dissociation processes by charged particles produced via three- and four-body 
decays of u. 

In order to derive typical constraint, in our analysis, the primary neutrino emitted from 
NLSP sneutrino is treated as equal- weight admixture of Ve-, i^fi, and Ur. Notice that the most 
important constraints, which are from hadro-dissociation processes, do not change even if 
we adopt different assumption (although the constraints from the two-body decay of i) may 
have slight dependence on the flavor of primary neutrino). Then, we compare the results of 
our calculation with light-element abundances inferred from observations; we calculate 
variable defined in |H] and obtained 95 % CL. constraints. The observational constraints 
adopted in this paper are summarized in Appendix 1X1 

Effects of long-lived sneutrino on BBN depends on primordial abundance of sneutrino. 
We parameterize the primordial abundance is parameterized by yield variable, which is 
defined as the ratio of number density and total entropy density (at t <^ To): 



With this quantity, the number density of sneutrino is given by n;^ = sY^e'^^'^'^ . Importantly, 
Yy depends on thermal history of the universe, so we consider several cases. 

We first consider the case where the sneutrinos are thermal relics; in this case, we assume 
"standard" evolution of the universe when the cosmic temperature is below the mass scale of 
SUSY particles msusY which is taken to be 100 GeV — 1 TeV. When the cosmic temperature 
T is higher than msusY? all the MSSM particles are thermalized and their number density 
is of order T^. As the temperature becomes lower than msusY? on the contrary, density of 
MSSM superparticles are Boltzmann-suppressed and the number densities of most of the 
superparticles become negligibly small as T — 0. Only the exception is the number density 
of the lightest superparticle in the MSSM sector, in our case, sneutrino. With i?-parity 
conservation, z/ decays only into gravitino and some other particle(s). Chemical equilibrium 
of sneutrino is maintained by pair annihilation process, whose rate becomes smaller than the 
expansion rate of the universe at some epoch. After this epoch, sneutrino freezes out from 
the thermal bath and its number density in the comoving volume is (almost) conserved at 
the cosmic time when t <^ tc,. Relic density of sneutrino in such a case was studied in [LlTj , 
and the yield variable is given by 



Using the yield variable given above, we calculate the abundances of light elements for fixed 
values of 777.3/2 and 777^. Then, comparing the results of our calculation with observational 
constraints on the primordial abundances of light elements, we derive constraints on 7713/2 
vs. rrip plane. 



Yu = 





(4.2) 
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Figure 4: Constrants on m^i2 vs. plane when the rehc sneutrinos have thermal origin. 
Regions inside the contours are excluded. 

In Fig.|3J we show constraints on 7713/2 vs. rriy plane when the relic sneutrinos have thermal 
origin. In this case, we have found that the most important constraints are obtained from 
the overproductions of D and ^Li. In order to understand the behavior of the constraints, it 
should be noticed that the lifetime becomes shorter as the gravitino mass becomes smaller. 
Consequently, when the gravitino mass is small enough, sneutrino decays at very early stage 
of BBN and no constraint is obtained. On the contrary, when 10^ sec ^Tp^ 10^ sec, the 
background "^He (which we denote obg) is effectively dissociated by the energetic hadrons 
produced in the hadronic shower. In this case, overproduction of D may occur as a result 
of hadro-dissociation of obg- In addition, energetic T and ^He is also produced by the 
hadro-dissociation process. Such T and ^He become sources of non-thermal production of 
^Li via T + obg ~^ ^Li + ^ and '^He + a^c ~^ ^Li + P- As the lifetime becomes longer, 
the energetic hadrons are stopped by the scattering processes with background particles and 
hence the effects of hadro-dissociations become inefficient. In addition, we also have found 
that effects of photo-dissociation are not important in this case. The constraint obtained 
here is qualitatively consistent with the result in ^Tj, but our constraint extends to the 
region with smaller gravitino mass. 

So far, we have considered the case where the relic sneutrinos have thermal origin; in 
such a case, the relic density of z/ is given by Eq. ()4.2|) . In non-standard cases, however, 
relic density of 9 may become larger. For example, for the case where late-decaying scalar 
condensation exists, superparticles may be produced at the decay time of scalar condensation. 
If the reheating temperature after the decay of scalar condensation is lower than ~ 10 GeV, 
relic density of sneutrino becomes larger than that given in Eq. ()4.2|) |18j . Thus, taking 
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Figure 5: Constraints on m^i2 vs. rriyYi) plane. Here, we take = 300 GeV. Dotted line 
shows the upper bound on m/jYp from the overproduction of gravitino due to the decay of 
sneutrino. 

as a free parameter, we calculate the abundances of light elements. 

In Fig. El we show constraints on 7713/2 vs. Yj> plane. Here, we take m,> = 300 GeV; 
with this value of sneutrino mass, the light-element abundances are consistent with obser- 
vational constraints if the relic sneutrinos have thermal origin. However, as we adopt larger 
value of Yy, some region of the parameter space is excluded. As one can see in Fig. El the 
most important constraints are from overproduction of ^He when 7713/2 ^0.4 GeV and from 
overproductions of D and ^Li when 7713/2 ^ 0.4 GeV. In addition, when the gravitino mass 
becomes close to the sneutrino mass, the lifetime of the sneutrino becomes relatively long and 
photo-dissociation of a-QG becomes effective. In this case, overproduction of ^He may also 
occur. On the other hand, for small gravitino mass 0.4 GeV) the lifetime of the sneutrino 
IS short (< 40 sec) and pion production by the high energy neutrino emitted in the two-body 
decay becomes significant, which leads to change p-n ratio and ^He overproduction. 

Gravitinos are produced by the decay of sneutrinos. Since gravitino is stable, they con- 
tribute to the present mass density of the universe. Density parameter of the gravitino 
produced by the sneutrino decay is obtained by 

^3/2 = "^"^'-^^ (4.3) 
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where Snow is the entropy density at the present universe and pc is the critical density. 
Requiring < ^cdm, is constrained from above. The bound is also shown in Fig. 
It should be noted that, when 777,3/2^0.! GeV, BBN constraints on are more stringent 
than that obtained from the overproduction of the gravitino in most of the cases. 



5 Implication for Leptogenesis 

In this section, we discuss implication of sneutrino-NSLP scenario for thermal leptogenesis. 
It is widely known that the present baryon asymmetry of the universe may originate from 
non-equilibrium decay of right-handed (s) neutrino which has frozen out from the thermal 

In order to realize thermal leptogenesis, however, high enough reheating temperature is 
needed j21j. Here, we define reheating temperature as 
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1/4 



T^^\ MtTfA , (5.1) 

\g*T^ j 

where is the effective number of massless degrees of freedom. (We take the MSSM value 

= 228.75 in our calculation.) Then, detailed calculation based on MSSM shows that 
Tr ^ 1 X 10^ GeV is required in order to generate large enough baryon asymmetry 

If the reheating temperature is as high as ~ 10^ GeV, however, large number of grav- 
itinos are produced by scattering processes of MSSM particles in thermal bath. Such grav- 
itinos usually cause serious cosmological problems. If gravitino is unstable, various colored 
and/or charged particles are produced when gravitinos decay, resulting in serious dissocia- 
tions of light elements during/after BBN. If Tr ~ 10^ GeV, light-element abundances are 
too much affected to be consistent with observations unless gravitino mass is extremely large 
(7713/2^10 TeV) |H101_. 

With stable gravitino, on the contrary, the present mass density of gravitino may become 
too large ^U]. Since the number density of gravitino is approximately proportional to the 
reheating temperature, we obtain upper bound on the reheating temperature. This fact 
severely constrains the scenario of thermal leptogenesis Importantly, as the gravitino 
mass becomes smaller, production of longitudinal mode of gravitino is enhanced and upper 
bound on Tr becomes more stringent. 

We have numerically solved Boltzmann equations from the infiaton-dominated epoch 
to radiation dominated epoch and derived upper bound on the reheating temperature. In 
our calculation, we have assumed that, in the infiaton-dominated epoch, infiaton potential 
is well approximated by parabolic potential and that infiaton is rapidly oscillating. Then, 



T^^In fact, gravitinos are also produced by scattering processes of particles in thermal bath. For details, 
see the next section. 

'^^Reheating temperature used in [^S] (which we denote j'^'^'^^^^)) jg related to our definition as 

y(GNRRS) ^ 
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Figure 6: Contours of constant maximal possible reheating temperature after inflation (or 
any other entropy production) in order to avoid overproduction of gravitino. The horizontal 
axis is the gravitino mass while the vertical axis is bino mass. 



denoting energy densities of radiation and inflaton as prad and pinf, respectively, and the 
number density of gravitino as ^3/2, they obey the following equations 

4ifprad + TinfPinf, (5.2) 
Sifpinf - TinfPinf, (5.3) 

3Hn3/2 + {atotVrc\)n^s,d, (5.4) 

where H is the expansion rate of the universe, and n^ad = ^i^^"^. In addition, (atot^rei) 
is "thermally averaged" total cross section (times relative velocity), which is given in P51. 
When the gravitino mass is small, (atot^rei) is proportional to squared of gaugino masses 
and, consequently, the resultant upper bound is sensitive to the gaugino masses. Here, we 
adopt GUT relation among gaugino masses 

Ms _ M2 _ 3 Ml 

9l ~ 9l~5 gl ' ^'-'^ 

where M3, M2, Mi and (73, (72, Qi are gaugino masses and gauge coupling constants for 
5'f/(3)c, SU{2)l, and U{1)y gauge interactions, respectively. 

Numerically solving Eqs. ()5.2|) — ()5.4|) from the epoch with H ^ Finf to the present 
epoch, we calculate the density parameter of gravitino ^3/2 = ^^^'"'^^^ as a function of 



"Prad 

dt 

dt 
dns/2 

dt 
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reheating temperature.*^ In deriving upper bound on Tr, we require that ^3/2 be smaller 
than the best-fit value of the dark- matter density reported by WMAP (see Eq. in 
other words, the upper bound T^^^^ satisfies the relation ^3/2(7^™^'''')/;,^ = 0.105. Thus, for 
the case where Tr = gravitino dark matter is realized. 

The result is shown in Fig. IHl As we mentioned, the gravitino production cross section 
is more enhanced with larger value of gaugino mass and smaller value of m3/2. Thus, in 
those limits, upper bound on Tr becomes more stringent. Comparing Fig. |3]with Fig. IHl 
we can see that reheating temperature of 0(10^ GeV) is possible without conflicting BBN 
constraints for 20 GeV ^777.3/2 ^400 GeV and Mi ^ 400 GeV. (Notice that m/) is smaller 
than Ml since consider the case where sneutrino is the lightest superparticles in the MSSM 
sector.) Thus in the scenario with gravitino NLSP and gravitino LSP, thermal leptogenesis 
scenario is possible in the above-mentioned parameter region. In this region, in addition, it 
should be also noticed that gravitino is a viable candidate for dark matter. 

It is important to note that the upper bound given in Fig. IHl is applicable irrespective 
of what the NLSP is. Thus, for thermal leptogenesis, gravitino mass should be 0(10 GeV) 
or larger. This fact gives serious constraint. When bino or charged slepton is the NLSP 
(with gravitino being the LSP), photo-dissociation processes provide stringent constraints 
on the parameter region 20 GeV < mg/a < 400 GeV and Mi < 400 GeV where Tr > 10^ GeV 
is allowed In particular, photo-dissociation processes of ^He enhance the ratio 

^He/D too much. Thus, if we adopt upper bound on the ratio ^He/D given in Eq. ()A.2|1 . 
such parameter region is excluded. Thus we have to conclude that thermal leptogenesis is 
difficult for the cases with bino NLSP and charged-slepton NLSP. On the contrary, as we 
have seen, thermal leptogenesis is possible in large parameter space if we consider sneutrino 
NLSP. Thus, the scenario with sneutrino NLSP and gravitino LSP has an advantage in 
realizing thermal leptogenesis. 

Before closing this section, we comment on what happens if we relax the GUT relation 
on the gaugino masses. With the GUT relation, gluino becomes much heavier than bino. 
Since the gravitino production cross sections are proportional to squared of gaugino masses, 
primordial gravitinos are mostly produced by processes with gluino in initial and/or final 
states. If the gluino mass is smaller, gravitino production rate is suppressed for a fixed value 
of Ml. Indeed, in some classes unified models, GUT relation may not hold; models with 
hypercolor is one of such cases j2Hl- Then, with smaller gluino mass, upper bound on 
the reheating temperature can become higher for a fixed value of bino mass. In such case, 
thermal leptogenesis may be possible even with bino NLSP and stau NLSP. 

6 Conclusions 

In this paper we have investigated the BBN constraint on the scenario in which the sneutrino 
is NLSP and decays into a gravitino (i.e., LSP) and the standard model particles. Sneutrino 

'^^In fact, gravitinos are also produce by the decay of sneutrino. However, we neglect such contribution 
since it is subdominant for most of the cases as as we can estimate from Eq. H4.2|) . 
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mainly decays into a gravitino and a neutrino which scatters off the background neutrinos 
and electrons producing pions and high energy electrons. Moreover quarks and photons 
are also produced via three- or four-body decay of the sneutrino although their branching 
ratios are small. For the case where the sneutrinos are thermal relics, the most stringent 
constraint comes from overproduction of D and ^Li produced in the hadron showers which 
are induced by the four-body decays. We also derived the constraint on the sneutrino 
abundance assuming the sneutrinos are produced non-thermally. We have found that the 
BBN constraint is more stringent than that from the cosmic density of the gravitino in the 
wide range of gravitino mass. It is also found that pion production by high energy neutrinos 
from the two-body decay becomes important when the sneutrinos decay at ^ 10 sec. 

Since the main decay mode {u ipv) leads to a less stringent constraint, the sneutrino 
NLSP scenario allows a larger parameter space than other scenarios like stau NLSP. In fact, 
it has been shown that the thermal leptogenesis is realized. 

Note Added: After finalizing the manuscript, we noticed the paper |^ which has some 
overlap with our analyses. 
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(K.K.), and the Grant-in-Aid for Scientific Research from the Ministry of Education, Science, 
Sports, and Culture of Japan, No. 15540247 (T.M.). 

A Observational Constraints on Light Elements 

In this appendix, we summarize the observational constraints on light-element abundances 
which we use in our analysis. For some light elements, we take account of recent developments 
in observations, so the observational constraints adopted in this paper is different from those 
used in the old study [71 |H] . 

The primordial value of the ratio D/H is measured in the high red-shift QSO absorption 
systems. Recently a new deuterium data was obtained from observation of the absorption 
system at the red-shift z = 2.70 towards QSO SDSS1558-0031 jHO]. The reported value 
of the deuterium abundance is log(D/H) = —4.48 ± 0.06. Combined with the previous 
data [SH E21 ISSl inn ISH] 5 it is reported that the primordial abundance is given by 



Here we have used the fact that the ^He/D is a increasing function of the time since D is 
more destroyed than ^He during chemical evolution, and hence the solar abundance gives 
the upper limit. 

The ^He is observed in metal poor extragalactic HII regions and the primordial abundance 
Yp is obtained by extrapolation of the observed abundances to zero metallicity. For a long 



For the constraint on 



(D/H)p = (2.82 ± 0.26) x 10^^ 
^He, we adopt the same constraint as in [H |H] : 
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time relatively low abundance had been believed for '^He. Field and Olive jB^ derived 
Yp = 0.238 ± (0.002),tat ± (0.005)syst and Izotov and Thuan ^ obtained Yp = 0.242 ± 0.002. 
However, recently, Olive and Skillman [SB] reanalyzed the Izotov-Thuan data and derived 
higher ^He abundance with much larger uncertainty, Yp = 0.249 ± 0.09. The similar value is 
obtained in more recent analysis jSHl which we adopt here, 

Fp = 0.250 ± 0.004. (A.3) 

The smaller uncertainty than that in [2Bj is due to different determination of the electron 
temperature Te for HIT regions. The OIII lines are used in jSH] while ^2SJ uses Hel recombi- 
nation lines, which generally leads to large errors in T^. 

Warmest metal-poor (pop. II) halo stars have almost constant ^Li abundances (Spite 
plateau) independent of metallicity. This constant value is considered as the primordial 
abundance of ^Li. Bonifacio et al. |4Q| obtained logio(^Li/H)p = —9.66 ±0.056 which is con- 
sistent with more recent value log]^o(^-'-'i/H)p = —9.63 ± 0.06 by Melendez and Ramirez |41] . 
On the other hand, significantly small abundance log^o(''Li/H)p = —9.91 ± 0.10 was derived 
in 02] where it was claimed that there is a correlation between ''Li and Fe abundances due 
to ^Li production by cosmic ray interactions and one should obtain the primordial value by 
taking zero Fe limit. The Fe dependence of ''Li was also observed by Asplund et al [121 who 
reported log]^Q(^Li/H)p = —9.90 ± 0.06. However, no correlation was observed in jlOl HT] . 
Thus, the ^Li-Fe correlation is still an open question. Furthermore, the observed ^Li abun- 
dance may be smaller than the primordial value if ''Li is depleted in stars. For example, 
03] shows that rotational mixing leads to depletion factor Dj at most 0.3 dex. Therefore, 
at present it is difficult for us to reach some consensus. Since the upperbound of ^Li is im- 
portant in deriving a constraint on the sneutrino decay, we conservatively adopt the higher 
value IH] and add systematic error of 0.3 dex taking depletion into account, 

logio(^Li/H)p = -9.63 ± 0.06 ± 0.3. (A.4) 

Here we have added —0.3 for systematic error because the high ^Li abundances in |HH l¥T] 
are different from those in |121l3ni by about factor of 2 (^ 0.3 dex). 

As for ^Li recently Asplund et al. |13] detected ^Li/^Li ratios in 9 metal poor halo dwarfs. 
In particular, ^Li abundance detected in very metal poor star LP 815-43 with [Fe/H] = - 
2.74 was is ^Li/''Li = 0.046 ± 0.022. Such high ^Li abundance seems difficult to explain 
by the Galactic cosmic ray spallation and a-fusion reactions and might be primordial. In 
this paper we regard it as upper bound on the primordial ^Li abundance. Moreover, the 
depletion in stars is more important for ^Li since ^Li is more fragile than ''Li. The depletion 
factor Dq for ^Li is related to Dj as Dq ~ 2.5Dt jl3], which leads to log]^Q(^Li/^Li)p = 
1.5Di + log^Q(^Li/^Li)obs Taking account of the depletion, we adopt 

(6Li/^Li)p < 0.046 ± 0.022 + 0.084. (A.5) 
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